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Abstract. Transforming constraint models is an important task in re- 
cent constraint programming systems. User-understandable models are 
defined during the modeling phase but rewriting or tuning them is manda- 
tory to get solving-eflicient models. We propose a new architecture al- 
lowing to define bridges between any (modeling or solver) languages and 
to implement model optimizations. This architecture follows a model- 
driven approach where the constraint modeling process is seen as a set 
of model transformations. Among others, an interesting feature is the def- 
inition of transformations as concept-oriented rules, i.e. based on types 
of model elements where the types are organized into a hierarchy called 
a metamodel. 



1 Introduction 

Constraint programming (CP) systems must combine a modeling language and 
a solving engine. The modeling language is used to represent problems with vari- 
ables, constraints, or statements. The solving engine computes assignments of 
variables satisfying the constraints by exploring and pruning the space of poten- 
tial solutions. This paper considers the constraint modeling process as constraint 
model transformations between arbitrary modeling or solver languages. It fol- 
lows several important consequences on the architecture of systems and user 
practices. 

Constraint programming languages are rich, combining common constraint 
domains, e.g. integer constraints or linear real constraints, with global constraints 
like alldif f erent, and even statements like if-then-else or forall. More- 
over the spectrum of syntaxes is large, ranging from computer programming 
languages like Java or Prolog to high-level languages intended to be more human- 
comprehensible. This may be contrasted with the existence of a standard lan- 
guage in the field of mathematical programming, which improves model sharing, 
writing and understanding. The quest of a standard CP language is a recent 
thread, dating back to the talk of Puget [15]. Another important concern is to 
employ the best solving technology for a given model. As a consequence, a new 
kind of architecture emerged. The key idea is to map models written with a 
high-level CP language to many solvers. For instance within the G12 project. 



MiniZinc [13] is intended to be a standard modeling language, and Cadmium [3] 
is able to map MiniZinc models to a set of solvers. Essence [5] is another CP 
platform offering an high level modeling language refining Essence specifications 
to Essence' models using Conjure [6] . Then hand- written translators can gener- 
ate models for several different solvers. The role of a mapping tool is to bridge 
modeling and solver languages and to optimize models for improving the solv- 
ing process. Cadmium is based on Constraint Handling Rules [8] and is the the 
closest CP platform from our model-driven approach. 

In our approach, we suppose that any CP language can be chosen at the 
modeling phase. In fact, finding a standard language is hard and existing lan- 
guages have their own features. It then becomes necessary to define mappings 
between any (pure modeling or solver) languages. This is just the first goal of 
the new architecture for constraint model transformations defined in the sequel. 
It follows many advantages: 

— Any user may choose its favourite modeling language and the known best 
solving technology for a given problem provided that the transformation 
between languages is implemented. 

~ It may be easy to create a collection of benchmarks for a given language 
from different source languages. This feature may speed up prototyping of 
one solver, avoiding hand rewriting of problems into the solver language. 

— A given problem may be handled using different solving technologies. Users 
may not have to play with solver languages. 

To this end, we define a generic and flexible pivot model (i.e. an intermediate 
model) to which any language is mapped. Considering a new language in this 
framework only requires a parser and a generally simple transformation to the 
pivot model. 

The second goal is to define refactoring operations and optimizations of con- 
straint models using declarative rules. Implementing them over pivot models 
guarantees the independence from external languages. In other words every op- 
eration is implemented once, by means of a so-called concept-oriented rule. In our 
model engineering approach the elements of models are specified within meta- 
models, which can be seen as a hierarchy of concepts or types. The rules are 
able to filter models according to these types, which may be more powerful than 
syntax-oriented rules. 

The third goal is to apply the best transformations for given solving tech- 
nologies. For instance, a matrix with a few non null elements could be trans- 
formed into a sparse matrix when using a linear algebra package. The selection of 
transformation steps is implemented as a sequential procedure, applying trans- 
formations until at least pivot models fit the structure requirements of the target 
language. 

This architecture has been fully implemented using a model-driven engineer- 
ing (MDE) approach [14]. MDE tools enable us to separate the grammar con- 
cerns from modeling concepts using dedicated tools and languages like TCS [11] 
and ATL [12, 10]. The main advantage is that we can reason about concepts and 
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Fig. 1. General transformation framework. 



their relations through a metamodel. Transformations are specified by defin- 
ing matchings between concepts at the metamodel level of abstraction. Thus, 
grammar concerns are relegated into the foreground, while concepts processing 
becomes the major task. 

With respect to previous works, e.g. [4], the new architecture gives more 
freedom in constraint modeling. s-COMMA is not always the source modeling 
language and refactoring steps can be chosen. Thus, users can play with any 
modeling language, until it is mapped to our platform. Dealing with a solver does 
not require to manipulate its language. Moreover, handling a new language or a 
new transformation in the system requires a few work. The main limitation of our 
approach is that only the modeling fragments of languages can be processed i.e., 
the declarative part. It is not possible to partially execute a computer program 
that builds the constraint store. 

This paper is organized as follows. Section 2 presents an overview of our 
general transformation framework. Next section introduces the metamodels of 
two CP languages illustrated on a well-known problem. The pivot metamodel 
and the transformation rules are introduced in Section 4. Section 5 presents 
the whole model-driven process including the possibility of selecting relevant 
mappings. The related work and a conclusion follow. 

2 The Model-Driven Transformation Framework 

Figure 1 depicts the architecture of our model-driven transformation framework, 
which is classically divided in two layers Ml and M2 [14]. Ml holds the models 
representing constraint problems and M2 defines the semantic of Ml through 
metamodels. Metamodels describe the concepts appearing in models, e.g. con- 
straint, variable, or domain, and the relations among these concepts, e.g. inher- 
itance, composition, or association. In this framework, transformation rules are 
defined to perform a complete translation in three main steps: translation from 
source model A to the pivot model, refactoring/optimization on the pivot model, 
and translation from the pivot model to target model B. Models A and B may 
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Fig. 2. Extract of the s-COMMA metamodel. 



be defined through any CP languages. The pivot model may be refined several 
times in order to adapt it to the desired target model (see Section 4). 

A main feature resulting from a model-driven engineering approach is that 
transformation rules operate on the metamodel concepts. For instance, unrolling 
a f orall loop is implemented once over the f orall concept, which is indepen- 
dent from the many syntaxes of f orall in CP languages. In fact, no grammar 
specification is required for the pivot model. Syntax specifications of CP lan- 
guages must be defined separately using specific tools achieving text-to-model 
or modcl-to-tcxt mappings like TCS [11], which implement both tasks. 



3 A Motivating Example 

In this section, we consider two CP languages, and we motivate the needs and 
the means for implementing transformations between them. 

ECL'PS'^ [17] is chosen as a leading constraint logic programming system. 
s-COMMA [16] is an object-oriented constraint language developed in our team. 
Their mctamodels are partially depicted in Figure 2 and 3 using UML class 
diagram notation. The roots of these hierarchies are equivalent, such that the 
model concept represents the complete constraint problem to be processed. 

In s-COMMA, a model is composed of a collection of model elements. A model 
element is either an enumeration, or a class, or a constant. Each class is com- 
posed of a set of class features which can be specialized in variables, constant or 
constraint zones. Variable with a type defined as a class is an object. Constraint 
zones are used to group constraints and other statements such as conditionals 
and loops. The concepts of global constraints and optimization objective are not 
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Fig. 3. Extract of the ECL'PS'= metamodel. 



shown here, but can be also defined. The concept of expressions are not detailed 
in this paper since it is based on classical operatored expressions using boolean, 
set and arithmetic operators. 

In the ECL'PS'^ metamodel, we propose to define a model as a collection 
of predicates holding predicate elements and variables. Predicate elements are 
variable features or statements. Variables features is either a constant value 
assignment, a domain definition, an array or a set definition related to a variable. 
In fact, we consider that variables arc implicitly declared through their features. 

Considering the well-known problem of the social golfers. Figure 4 and 5 show 
two versions of the same problem using s-COMMA and ECL'PS'^ languages. This 
problem considers a group of n = g x s golfers that wish to play golf each week, 
arranged into g groups of s golfers, the problem is to find a playing schedule for 
w weeks such that no two golfers play together more than once. 

The s-COMMA model is divided in a data file and a model file. The data file 
contains the golfer names encoded as an Enum concept at line 1 and the problem 
dimensions defined by means of constants (size of groups, number of weeks, and 
groups per week). The model file represents the generic social golfers problem 
using the Model concept. The problem structure is captured by the three classes 
SocialGolf ers. Group, and Week, which are conformed to the Class concept. The 
Group class owns the players attribute corresponding to a set of golfers playing 
together, each golfer being identified by a name given in the enumeration from 
the data file. In this class, the constraint zone groupSize (lines 30 to 32) restricts 
the size of the golfers group. The Week class has an array of Group objects and the 
constraint zone playOncePerWeek ensures that each golfer takes part of a unique 
group per week. Finally, the SocialGolf ers class has an array of Week objects 
and the constraint zone dif f erentGroups states that each golfer never plays two 
times with the same golfer throughout the considered weeks. 
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// Data file 


15 


} 


2 


enum Name :— {a,b,c,d,e,f, 


16 


} 




g ,h, i }; 


17 


} 


3 


int s :— 3; // size of 


18 


class Week { 




groups 


19 


Group groups [ g ] ; 


4 


int w :— 4; //number of 


20 


constraint 




weeks 




pi ay Once Per Week { 


5 


int g :— 3; //groups per 


21 


forall{gl in l..g) { 




week 


22 


forall(g2 in gl + l..g) { 


6 




23 


card { groups [gl]. 
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// Model file 




players intersect 


8 


main class SocialGolfers { 




groups [g2]. players 


9 


Week weeks [w] ; 




) = 0; 


10 


constraint 


24 


} 




differentGroups { 


25 


} 


11 


forall(wl in l..w) { 


26 


} 


12 


forall(w2 in wl + l..w) { 


27 


} 


13 


card ( weeks [wl ] . groups [ 


28 


class Group { 




gl ] . players 


29 


Name set players; 




intersect weeks [w2 


30 


constraint groupSizc { 




] . groups [ g2 ] . 


31 


card (players) — s; 




players )<— 1 ; 


32 


} 


14 


} 


33 


} 



Fig. 4. The social golfers problem expressed in s-COMMA. 
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( for (GROUPS, 1 ,G) , param( 
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i n t s e t s ( 
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L = WEEKS.GROUPSJLAYERS, 


24 


nth(V7,V6.L) , 
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#(V7, V8) , V8 $= S , 
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(for(Wl,l ,W) , param(L,W,G 
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) do 
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(for(Gl,l ,G) ,param(L,G, 
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( for (W2,W1+1,W) , param(L 
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( for (G2, Gl + l.G) .param( 


10 


(for(Gl,l ,G) , param(L,G 




L,G, WEEKS, Gl) do 




,W1,W2) do 
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V9 is G* (WEEKS- 1)+G1, 
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( for (G2, 1 ,G) , param(L, 
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nth(V10,V9,L) , 




G,W1,W2,G1) do 


31 


Vll is G* (WEEKS- 1)+G2 


12 


VI is G*(W1-1)+G1, 






13 


nth(V2,Vl,L) , 


32 


nth(V12.Vll,L) , 


14 


V3 is G*(W2-1)+G2, 


33 


#(V10 /\ V12, 0) 


15 


nth(V4,V3,L) , 


34 
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#(V2 /\ V4, V5) ,V5 
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) 




$=< 1 
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39 


label_scts (L) . 



Fig. 5. The social golfers problem expressed in ECL'PS'^. 



Figure 5 depicts the ECL^PS"^ model resulting from an automatic transfor- 
mation of the previous s-COMMA model. The problem is now encoded as a single 
predicate whose body is a sequence of atoms. The sequence is made of the prob- 
lem dimensions, the list of constrained variables L, and three statements resulting 
from the transformation of the three s-COMMA classes. It turns out that parts 
of both models are similar. This is due to the sharing of concepts in the un- 
derlying mctamodels, for instance constants, f orall statements, or constraints. 
However, the syntaxes are different and specific processing may be required. For 
instance, the f orall statement of ECL*PS*^ needs the param keyword to declare 
parameters defined outside of the current scope, e.g. the number of groups G. 

The treatment of objects is more subtle since they must not participate to 
ECL^FS"^ models. Many mapping strategies may be devised, for instance map- 
ping objects to predicates [16]. Another mapping strategy is used here, which 
consists in removing the object-based problem structure. Flattening the prob- 
lem requires visiting the many classes through their inheritance and composi- 
tion relations. A few problems to be handled are described as follows. Impor- 
tant changes on the attributes may be noticed. For example, the weeks array of 
Week objects defined at line 9 in Figure 4 is refactored and transformed to the 
WEEKS_GRQUPS_PLAYERS fiat list stated at line 5 in Figure 5. It may be possible to 
insert new loops in order to traverse arrays of objects and to post the whole set 
of constraints. For instance, the last block of for loops in the ECL'PS'^ model 
(lines 27 to 39) has been built from the playOncePerWeek constraint zone of the 
s-COMMA model, but there is two additional for loops (lines 21 and 22) since the 
Week instances arc contained in the weeks array. Another issue is related to lists 
that cannot be accessed in the same way than arrays in s-COMMA. Thus, local 
variables (v^) and the well-known nth Prolog built-in function are introduced in 
the ECL^PS-^ model. 

4 Pivot metamodel and refactoring rules 

The pivot model of a constraint problem is an intermediate model to be trans- 
formed by rules. The rules may be chained to implement complex transforma- 
tions. In the following, the pivot and some structural refactoring and optimiza- 
tion rules are presented. 

4.1 Pivot metamodel 

Our pivot model has been designed to support as much as possible the constructs 
present in CP languages, for instance variables of many types, data structures 
such as arrays, record, classes, first-order constraints, common global constraints, 
and control statements. We believe that it is better and simpler to establish a 
general CP metamodel, while it is more complex to find a standard CP concrete 
syntax. 

Figure 6 depicts the metamodel associated to pivot models. A pivot model 
is composed of a collection of elements, divided in three main concepts: types. 
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Fig. 6. Extract of the pivot metamodel. 



features and the concrete concept of predicate. The inheritance tree of types is 
the same as in the s-COMMA metamodel (see Figure 2). The inheritance tree for 
model features is also quite similar, except for the concept of record which is an 
untyped collection of features. 



4.2 Pivot model refactoring 

We define several refactoring steps on pivot models in order to reduce the pos- 
sible gap between source and target model. These steps are implemented in 
several model transformations, most of them being independent from the oth- 
ers. The idea is to refine and optimize models in order to fit the target languages 
supported concepts. 

Model transformations are implemented in the declarative transformation 
rule language ATL [12]. This rule language is based on a typed description of 
models to be processed, namely their metamodel. In this way, rules are able to 
clearly state how concepts from source metamodels are mapped to concepts from 
the target ones. For the sake of simplicity, only a few of the more representative 
rules of transformations are shown. ATL helpers are not detailed, but they only 
consist of OCL navigation. 



Composition flattening This refactoring step replaces object variables by 
duplicating elements defined in their class definition. Names of duplicated vari- 
ables are prefixed using their container name in order to avoid naming ambi- 
guities. This refactoring step processes object variables and their occurrences, 
while other entities are copied without modification. In fact, two ATL transfor- 
mations are defined to ease each refactoring step. The first one removes classes 
and object variables by replacing them by the concept of record (see Figure 7). 



It can be highlighted that there is no ATL rule where the source pattern matches 
elements being instances of CSPClass. Thus, they are implicitly removed from 
models (obviously no rule creates class instances). The second transformation 
removes records to get flattened variables (see Figure 8). 



1 rule Model { 

2 from 

3 s : Pivot ! CSPModcl 

4 to 

5 t : Pivot ! CSPModcl ( 

6 name <— s . name , 

7 elements <— s . elements 

8 } 

9 rule Variable { 

10 from 

11 s : PivotCSP! CSPVariable ( 

12 not s . inustBeDuplicated 

13 ) 

14 to 

15 t : PivotCSP! CSPVariable ( 

16 name <— s . name , 

17 type <— s . type , 

18 domain <— s . domain , 

19 i s S e t <— s . i s S e t , 

20 array <— s . array 

21 ) 

22 } 

23 rule Variablc2Rccord { 

24 from 

25 s : PivotCSP! CSPVariable ( 

26 s.isObject 

27 ) 

28 to 

29 t : PivotCSP! CSPRccord ( 

30 name <— s . name , 

31 array <— s . array , 

32 elements <— s . type . features — > collect (f| 

33 thisModule.duplicate(f) 

34 ) 

35 ) 

36 } 



Fig. 7. An extract of ATL rules used to remove the concept of class in pivot models. 



In Figure 7, the first rule (lines 1 to 8) is used to copy the root concept of 
model. Most of other concepts are duplicated with similar rules like the the 
second one (lines 9 to 22). The helper mustBeDuplicated is defined for each 
CSPModelFeature and it returns true when: (1) the considered element is an ob- 
ject variable (its type is a class) or (2) it is a feature of a class. Using the last 
rule, object variables are replaced by records. The helper isObject returns true 
only if the type of variables is a class. In this rule, features of variable classes 
are browsed using OCL navigation (collect statement over s. type. features). 
The rule duplicate is applied on each feature. This rule is lazy and abstract. It 
is specialized for each CSPModelFeature concrete sub-concepts and it creates as 
many features as it is called. 



The second transformation processes records by replacing them by their set 
of elements. This is easily done by collecting their elements from their container 
as shown on Figure 8 at lines 7 to 11. The helper getAllElements returns the set 
of CSPModelFeature within a record or a hierarchy of records. 



1 rule CSPModel { 

2 from 

3 s : PivotCSPICSPModol 

4 to 

5 t : PivotCSPICSPModol ( 

6 name <— s . name , 

7 elements <— s . elements — > union ( s . elements — > select (r j 

8 r . ocllsTypeOf (PivotCSP ! CSPRocord) 

9 )->collect (r | 

10 r. getAllElements 

11 )->flatten () ) 

12 ) 

13 } 

14 rule RccordArray { 

15 from 

16 s : PivotCSP! CSPRecord ( 

17 (not s . ar r ay . o cl 1 s U n d e f i n cd ( ) ) and 

18 s. elements — > select (c| 

19 e . ocllsKindOf(PivotCSP! CSPStatcmcnt ) 

20 )->size()>0 

21 ) 

22 to 

23 t : PivotCSP! CSPForall ( 

24 index <— i , 

25 constraints <— s . elements — > reject (c 

26 e . oclIsKindOf(PivotCSP ! CSPTypedElement) 

27 ) , 

28 i : PivotCSP! CSPlndexVariable ( 

29 name <— s . name , 

30 domain <— d 

31 ) , 

32 d : PivotCSP! CSPIntcrvalDomain( 

33 lower <— 1 , 

34 upper <— thisModule . duplicateExpr( s . array . n) 

35 ) , 

36 1 : PivotCSP! CSPIntVal( 

37 value <- 1 

38 ) 

39 } 



Fig. 8. Main ATL rules used to remove the concept of record in pivot models. 



However, some other complex rules must be defined to process arrays of 
records, (formerly arrays of object variables). Indeed, contained statements have 
to be encapsulated in a for loop to take into account the constraints for all objects 
in the array. This task is performed by the rule RecordArray which create a new 
for loop over the record statements (lines 25 to 27). A new for loop requires also 
a new index variables with its domain (lines 28 to 38). 

Using the concrete syntax of s-COMMA, Figure 9 shows the result of this 
refactoring step. The name of the variable at line 1 corresponds to the concate- 



nation of all object variable names. The two for loops (lines 2 and 3) were created 
from the arrays of objects using their name for index variables. 
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int set wccks_groups_playcrs [ w* g ] in [ 1 , 9 ] , 
forall weeks in [1 ,w] { 
forall groups in [l,g] { 

card (weeks_groups_players [weeks *w-f-gr oup s ] ) — g , 



} 



} 



Fig. 9. Extract of the social golfers pivot model after composition removal and enu- 
meration removal transformations. 



Enumeration removal During this refactoring step, enumeration variables are 
replaced by integer variables with a domain defined as an interval from one to 
the number of elements within the enumeration. Line 1 in Figure 9 shows the 
result of this transformation on the enumeration called Name in the social golfers 
model: the variable has an integer domain from 1 to 9 replacing the set of nine 
values {a,b,c,d,e,f,g,h,i}. In the same way, occurences of CSPEnumLiteral are 
replaced by their position in the sequence of elements of the enumeration type. 

Other implemented refactoring steps Some other generic refactoring steps 
have been implemented in ATL to handle some structural needs. They arc not 
detailed since their complexity is similar to the previous examples and to detail 
all of them is not the scope of this paper. 

— If statements can be replaced by one constraint based on one or two boolean 
implications. For instance, if a then b else c becomes (a — 6) A (-la — > c). 

— Loop structures can be unrolled, i.e. the loop is replaced by the whole set of 
constraints it implicitly contains. Within expressions, the iterator variable 
used by the loop structure is replaced by an integer corresponding to the 
current number of loop turns. 

— Expressions can be simplified if they are constants. Boolean and integer 
expressions are replaced by their evaluation. Real expressions are not pro- 
cessed, because of real number rounding errors. More subtle simplifactions 
can be performed on boolean expressions such as a V -la that is always true. 
Only atomic boolean elements are processed by this last step. 

— Matrices are not allowed in all CP language, thus they can be replaced by one 
dimension arrays. Their occurrences in expressions must also be adapted: the 
index of the array is computed as follows: j] becomes m[j + (i* ncols)], 
where ncols is the number of columns of the matrix to. 

— The ECL'PS*^ language does not allow some sort of expressions. For instance, 
arrays of int sets cannot be accessed like other arrays with ' [ ] ' . Thus, an 



ECL'PS"^ specific transformation processes expressions and introduces focal 
variabfes if needed, as siiown on Figure 5 witii variabfes and nth predicate 
caffs. 

5 Handling CP languages and transformation chains 

In tfiis section, we describe the whofe transformation cfiain from a given CP 
fanguage to anotfier fanguage. 

5.1 Parsing CP languages 

The front-end of our system parses a source CP fanguage fife to get a modef rep- 
resentation (on which transformation rufes act) matching the concepts of the CP 
fanguage (injection phase). The back-end generates the code in the target CP 
fanguage (extraction phase) from the modef representation. Interfacing CP fan- 
guages and metamodefs is impfemented by means of the TCS toof [11]. This tool 
affows one to smoothfy associate grammars and metamodefs. It is responsibfe 
for generating parsers of CP fanguages and afso code generators. 

Figure 10 depicts an extract of the TCS fife for s-COMMA. In a TCS fife 
every concrete concept must have a corresponding template to be matched. For 
instance, the SCMAClass tempfate impfements the grammar pattern for cfass dec- 
farations using at the same time features of this concept defined in the meta- 
modef of s-COMMA. At parsing time on tfie s-COMMA sociaf goffers exampfe 
(see Figure 4, the "class" tolcen is matched for the week; cfass statement. Then 
Week is processed as the name attribute (a string in tfic metamodef) of a new 
class instance. Then the "{" token is recognized and the class features (the ar- 
ray of groups and the constraint) are processed by impficit matchings to their 
corresponding templates using the features reference. Finally the "}" token ter- 
minates the pattern description. In the SCMAClass template (lines 4 to 8), several 
TCS keywords are used. Here is a description of the most important keywords 
use in Figure 10: 

— context defines a local symbol table. 

— addToContext adds instances to the current symbol table. 

— refersTo accesses to the symbol tables according to the given parameter (here 
the name) to check the existence of an already declared element. 



5.2 Model checking rules 

The presented mctamodcls (see section 2) and the previous subsection show how 
to get CP language models. However, many irrelevant or erroneous models can 
be obtained without any additional checking [2]. For instance, variables may 
be defined with empty domains or expressions may be ill made (e.g. several 
equalities in an equality constraint). 
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template SCMAModcl main context 

: elements; 
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9 
10 
11 
12 
13 
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template SCMAClass context addToContext 

: (isMain ? "abstract main'') (isAbs tract ? "abstract") 
class" name 

(isDefined(supcrTypes) ? "extends" 

superTy pes { separator —" ," , refersTo— name , 
importContext } ) 
"{" [ features {separator—"."} ] "}" ; 



template SCMAVariable addToContext 

: (isSet ? "set" : "") type { refersTo— name} 
name ( isDefined ( array ) ? array) 

( isDefined { domain ) ? "in" domain); 



Fig. 10. Linking the grammar and the metamodel of s-COMMA with TCS. 

Several ATL transformations are used to check source models. We transform 
a source CP model to a model conform to the metamodel Problem defined in 
the ATL zoo'^. A Problem model corresponds to a set of Problem elements. This 
concept is only composed of three features: 

— severity is an attribute with an enumerated type which possible values are: 
error, warning and critic. 

— location is a string used to store le location of the problem in the source 



— description is a string used to defined a relevant message to descibe the 



Multiple ATL rules have been implemented to check models. Here is an ex- 
tract of the list of properties to check: 

— Some type checking on expressions. Operands must have a consistent type 
with the operator. For instance, an equality operator may operate on arith- 
metic expressions. 

— The consistency of variable domains : they must be based on constant ex- 
pressions and interval domains must have a lower bound smaller than the 
upper bound. 

— No composition or inheritance loops in s-COMMA. 

5.3 Chaining model transformations 

After the injection step or before the extraction step, models have to be trans- 
formed with respect to our pivot metamodel. All the refactoring steps presented 
in Section 4.2 are clearly not necessary in a transformation chain. Indeed, it 
clearly depends on the modeling structures of the source and target CP lan- 
guages. The idea is to use most of constructs supported by the target language 
to have a target model close, in terms of constructs, to our source model. For 



file. 



problem. 



^ http://www.eclipse.Org/m2m/atl/atlTransformations/#KM32Problem 



instance, when translating a s-COMMA model to ECL^PS"^, we should transform 
the objects. So, we choose the composition flattening step. We also need the 
enumeration removal and other refactoring steps such as the use of local vari- 
ables and nth predicates. Optionally, we may select the expression simplification 
steps. 

The whole transformation chain is based on three kind of tasks: (1) injec- 
tion/extraction steps, (2) transformation steps from/to the pivot model, (3) 
relevant refactoring steps. Transformation chains arc currently performed us- 
ing Ant scripts^. These scripts are hand- written, but they can be automatically 
generated using the am3 tool [1] and the concept of megamodel [7] to get a 
graphical interfaces to manage terminal models, metamodels and complex trans- 
formation chains. However, Automating the building of transformation chains is 
not possible with current tools. It would require to deeply analyze models and 
transformations to build relevant transformation chains. 

6 Experiments 

The benchmarking study was performed on a 2.66Ghz computer with 2GB RAM 
running Ubuntu. The ATL regular VM is used for all model-to-model transfor- 
mations, whereas TCS achieve the text-to-model and model-to-text tasks. Five 
CP problems were used to validate our approach as shown in Table 1. The sec- 
ond column represents the number of lines of the s-COMMA source files. The next 
columns correspond to the time of atomic steps (in seconds) : model injection (In- 
ject), transformations from s-COMMA to Pivot (s-to-P), refactoring composition 
structures (Comp), refactoring enumeration structures (Enum), transformations 
from Pivot to ECL'PS'^ (P-to-E), and target file extraction (Extract). The next 
column details the total time of complete transformation chains, and the last 
column corresponds to the number of lines of the generated ECL'PS'^ files. 



Problems 


Lines 


Inject 


s-to-P 


Comp 


Enum 


P-to-E 


Extract 


Total 


Lines 




(-) 


(s) 














(-) 


SocialGolfers 


42 


0.107 


0.169 


0.340 


0.080 


0.025 


0.050 


0.771 


38 


Engine 


112 


0.106 


0.186 


0.641 


0.146 


0.031 


0.056 


1.166 


78 


Send 


16 


0.129 


0.160 


0.273 




0.021 


0.068 


0.651 


21 


StablcMarriagc 


46 


0.128 


0.202 


0.469 


0.085 


0.027 


0.040 


0.951 


26 


10- Queens 


14 


0.132 


0.147 


0.252 




0.017 


0.016 


0.564 


12 



Table 1. Times for complete transformation chains of several classical problems. 



The transformation chain is efficient for these small problems. The text file 
injection and extraction are fast. The parsing phase is more expensive than the 
extraction, since it requires the management of symbol tables. The extraction 
phase settle for reading the ECL'PS'^ model. It can also be noticed that model 
transformations to and from the pivot are quite efficient, more especially the 

* http: / /wiki. eclipse.org/index.php / AM3_Ant_Tasks 



transformation to ECL'PS'^ model. It can be explained by the refactoring phases 
on the pivot model which simplify and reduce the data to process. We see that 
the composition flattening step is the more expensive. In particular, the Engine 
problem exhibits the slowest running time, since it corresponds to the design of 
an engine with more object compositions. 



Problems 


Inject 


s-to-P 


Comp 


Forall 


P-to-E 


Extract 


Total 


Lines 


Total/Lines 












(s) 




i^) 




(-) 


(-) 


5-Quccns 


0.132 


0.147 


0.252 


0.503 


0.071 


0.019 


1.124 


80 


«0.014 


10-C 


^uccns 


0.132 


0.147 


0.252 


1.576 


0.280 


0.060 


2.447 


305 


«0.008 


15-C 


Queens 


0.132 


0.147 


0.252 


3.404 


0.659 


0.110 


4.704 


680 


«0.007 


20-C 


Queens 


0.132 


0.147 


0.252 


6.274 


1.224 


0.178 


8.207 


1205 


«0.006 


50-C 


Queens 


0.132 


0.147 


0.252 


32.815 


13.712 


1.108 


48.166 


7505 


«0.006 


75-C 


Queens 


0.132 


0.147 


0.252 


80.504 


54.286 


2.456 


137.777 


16880 


«0.008 


100- 


Queens 


0.132 


0.147 


0.252 


175.487 


126.607 


4.404 


307.029 


30005 


wO.OlO 



Table 2. Time of complete transformation chains of the N-Queens problem. 



Table 2 presents seven different sizes of the N- Queens problem where the 
loop unrolling step has been applied. This experiment allows us to check the 
scalability of our approach according to model sizes. It can be analyzed through 
the ratio given in the last column which aims at quantifying the efficiency of a 
transformation chain considering the execution time per generated lines. 

As shown on this table, the ratio first decreases, but after 50-Queens it slowly 
grows up. In fact, the first four row ratios are impacted by the steps before 
the loop unrolling process, but for the last three rows they become neglectible 
comparing to the whole execution time. It may be noticed that for big problems 
(after 50-Queens) the ratio smoothly increases. We can thus conclude that our 
approach is applicable even for huge models, although translations times are not 
the major concerns in CP. 

7 Conclusion and Future Work 

In this paper, we propose a new framework for constraint model transformations. 
This framework is supported by a set of MDE tools that allow an easy design 
of translators to be used in the whole transformation chain. This chain is com- 
posed by three main steps: from the source to the pivot model, refining of the 
pivot model and from the pivot model to the target. The hard transformation 
work (refactoring/optimization) is always performed by the pivot which provide 
reusable and flexible transformations. The transformations from/to pivot be- 
come simple, thus facilitating the integration of new language transformations. 
In this paper, only two languages are presented, but translation processes with 
Gecode and Rcalpaver [9] are already implemented. 

In a near future, we intend to increase the number of CP languages our ap- 
proach supports. We also want to define more pivot refactoring transformations 
to optimize and restructure models. Another major outline for future work is to 



improve the management of complex CP models transformation chains. Models 
can be qualified to determine their level of structure and to automatically choose 
the required refactoring steps according to the target language. 
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Preface 



This textbook is intended for use by students of physics, physical chemistry, 
and theoretical chemistry. The reader is presumed to have a basic knowledge 
of atomic and quantum physics at the level provided, for example, by the first 
few chapters in our book The Physics of Atoms and Quanta. The student of 
physics will find here material which should be included in the basic education 
of every physicist. This book should furthermore allow students to acquire an 
appreciation of the breadth and variety within the field of molecular physics and 
its future as a fascinating area of research. 

For the student of chemistry, the concepts introduced in this book will provide 
a theoretical framework for that entire field of study. With the help of these con- 
cepts, it is at least in principle possible to reduce the enormous body of empirical 
chemical knowledge to a few basic principles: those of quantum mechanics. In 
addition, modern physical methods whose fundamentals are introduced here are 
becoming increasingly important in chemistry and now represent indispensable 
tools for the chemist. As examples, we might mention the structural analysis of 
complex organic compounds, spectroscopic investigation of very rapid reaction 
processes or, as a practical application, the remote detection of pollutants in the 
air. 
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Abstract. Tlie abstract sliould summarize the contents of the paper 
using at least 70 and at most 150 words. It will be set in 9-point font 
size and bo inset 1.0 cm from the right and left margins. There will be 
two blank lines before and after the Abstract. . . . 

1 Fixed-Period Problems: The Sublinear Case 

With this chapter, the prehminaries are over, and we begin the search for periodic 
solutions to Hamiltonian systems. All this will be done in the convex case; that 
is, we shall study the boundary-value problem 



with H{t, •) a convex function of x, going to +00 when — > 00. 
1.1 Autonomous Systems 

In this section, we will consider the case when the Hamiltonian H{x) is au- 
tonomous. For the sake of simplicity, we shall also assume that it is C^. 

Wc shall first consider the question of nontriviality, within the general frame- 
work of (Aoo, i?oo)-subquadratic Hamiltonians. In the second subsection, we shall 
look into the special case when H is (0, 6oo)-subquadratic, and we shall try to 
derive additional information. 

The General Case: Nontriviality. We assume that H is (Aqo, .Boo)-sub- 
quadratic at infinity, for some constant symmetric matrices A^o and B^o, with 
Boo — A-oo positive definite. Set: 



X = JH'{t,x) 
x{0) = x{T) 




7 : = smallest eigenvalue of Boo — A 



00 



(1) 



A : = largest negative eigenvalue of J— + A, 



dt 



■00 



(2) 



Theorem 1 tells us that if A + 7 < 0, the boundary-value problem: 

X = JH'{x) 
x(0) = x{T) 



(3) 



has at least one solution x, which is found by minimizing the dual action func- 
tional: 



{A-\,u)+N*{-u) 



dt 



(4) 



on the range of A, which is a subspace R{A)j^ with finite codimension. Here 



N{x) := H{x)--{A^x,x) 

is a convex function, and 

N{x) <^{{B^- Aoc)x,x)+c yx 
Proposition 1. Assume H'{0) = and H{0) = 0. Set: 

5 :=liminf2iV(a;) ||x|r^ . 

Vl < ~^ < ^! the solution u is non-zero: 

x{t) ^0 Vt . 



(5) 



(6) 



(7) 



(8) 



Proof. Condition (7) means that, for every 6' > 6, there is some e > such that 

(9) 



\\x\\<e^N{x)<^-\\xf 



It is an exercise in convex analysis, into which we shall not go, to show that 
this implies that there is an 77 > such that 



f\\x\\<ri^N*{y)<^\\y\f 



(10) 



Fig. 1. This is the caption of the figure displaying a white eagle and a white horse on 
a snow field 
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Since ui is a smooth function, we will have <ri for h small enough, 

and inequality (10) will hold, yielding thereby: 



i'ihui) < ~ \\ui\\l + \\uif 



(11) 



If we choose d' close enough to 5, the quantity + ^ill t)e negative, and 
we end up with 

< for ft 7^ small . (12) 

On the other hand, we check directly that tp{0) = 0. This shows that cannot 
be a minimizer of tjj, not even a local one. So u ^ and u ^ yl~^(0) =0. □ 

Corollary 1. Assume H is and {aoo,boo)-subquadratic at infinity. Let ^i, 
...,^7V be the equilibria, that is, the solutions of H'{^) = 0. Denote by cok the 
smallest eigenvalue of H" (^fe), and set: 



CO := Min {coi, . . . , cOk} 



If: 



T 



(13) 
(14) 



then minimization ofijj yields a non-constant T -periodic solution x. 



We recall once more that by the integer part E[a] of a € IR, we mean the 
a E Z such that a < q < a + 1. For instance, if we take aoo = 0, Corollary 2 
tells us that x exists and is non-constant provided that: 



or 



27r 277 

w ' boo 



(15) 
(16) 



Proof. The spectrum of yl is + aoo- The largest negative eigenvalue A is 

given by -\- floo, where 



27r, „ ^ 27r 

— fco + aoo < < — (fco + 1) + 



Hence: 



ko = E 

The condition 7 < —A < 6 now becomes 



aoo < 



27r 



fco - Ooo < w - ao 



which is precisely condition (14). 



(17) 
(18) 

(19) 

□ 
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Lemma 1. Assume that H is C2 on IR^"\{0} and that H"{x) is non-degenerate 
for any x ^0. Then any local minimizer xofip has minimal period T. 

Proof. We know that x, or x + ^ for some constant ^ G IR^", is a T-periodic 
solution of the Hamiltonian system: 

X = JH'{x) . (20) 

There is no loss of generality in taking ^ = 0. So ip{x) > ipix) for all x in 
some neighbourhood of x in VK^^^ (]R/T^; IR^") _ 

But this index is precisely the index ixix) of the T-periodic solution x over 
the interval (0,T), as defined in Sect. 2.6. So 

irix) = . (21) 

Now if X has a lower period, T/k say, wc would have, by Corollary 31: 

irix) = ihT/k{x) > kiT/k{x) + k-l>k-l>l . (22) 

This would contradict (21), and thus cannot happen. □ 

Notes and Comments. The results in this section are a refined version of [1]; the 
minimality result of Proposition 14 was the first of its kind. 

To understand the nontriviality conditions, such as the one in formula (16), 
one may think of a one-parameter family xt, T G (27ra;~^, 27r6^^) of periodic 
solutions, xt{0) = xt{T), with xt going away to infinity when T — >■ 2Tru)~^, 
which is the period of the linearized system at 0. 



Table 1. This is the example table taken out of The TgXbook, p. 246 



Year 


World population 


8000 B.C. 


5,000,000 


50 A.D. 


200,000,000 


1650 A.D. 


500,000,000 


1945 A.D. 


2. .300, 000.000 


1981) A.D. 


l,il)l),()()()AKK) 



Theorem 1 (Ghoussoub-Preiss). Assume H(t,x) is {0,e)-subquadratic at 
infinity for all £ > 0, and T-periodic in t 

H{t, ■) is convex Vt (23) 

H{-,x) is T-periodic Va; (24) 
H{t,x) > n{\\x\\) with n(s)s~"^ — >■ oo as s — >■ oo (25) 
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Ve>0, 3c : H{t,x) < ^\\xf + c . (26) 

Assume also that H is , and H"(t,x) is positive definite everywhere. Then 
there is a sequence Xk, fc € IN, of kT -periodic solutions of the system 

x = JH'{t,x) (27) 

such that, for every A; G IN, there is some Po G IN with: 

P>Po^Xpkj^Xk . (28) 

□ 

Example 1 (External forcing^. Consider the system: 

X = JH'{x) + f{t) (29) 

where the Hamiltonian H is (0, 6(x))-subquadratic, and the forcing term is a 
distribution on the circle: 

f=j^F + fo with e {M/TZZ- M^") , (30) 
where fo := jj f{t)dt. For instance, 

m = E '^'^^ ' (31) 

feeiN 

where Sk is the Dirac mass at t = fc and ^ G K^" is a constant, fits the pre- 
scription. This means that the system x = JH'(x) is being excited by a series 
of identical shocks at interval T. 

Definition 1. Let Aoo{t) and Boo{t) he symmetric operators m H^", depending 
continuously on t ^ [0,7"], such that Aoo{t) < Boo(^) for all t. 

A Borelian function H : [0,T] x H^" IR is called {A^,B^)-subquadratic 



at infinity if there exists a function N{t, x) such that: 

Hit, a;) = i {A^{t)x, x) + N{t, x) (32) 

Vt , N{t, x) is convex with respect to x (33) 

N{t,x) > n{\\x\\) with n{s)s~'^ ^ +00 as s +00 (34) 

3c G R : H{t, x)<^ iB^{t)x, x) + c \/x . (35) 



If Aao{t) = Gool and B^{t) = boci, with a^o < 600 G IR? we shall say that 

H is (aQc,boo)-subquadratic at infinity. As an example, the function with 
1 < a < 2, is {0, e)-subquadratic at infinity for every e > 0. Similarly, the 
Hamiltonian 

H{t,x) = ^k\\kf + \\xr (36) 
is {k, k + e)-subquadratic for every £ > 0. Note that, if k <0, it is not convex. 
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Notes and Comments. The first results on subharmonics were obtained by Ra- 
binowitz in [5] , who showed the existence of infinitely many subharmonics both 
in the subquadratic and superquadratic case, with suitable growth conditions 
on H'. Again the duality approach enabled Clarke and Ekeland in [2] to treat 
the; same problem in the convex-subquadratic case, with growth conditions on 
H only. 

Recently, Michalek and Tarantello (see [3] and [4] ) have obtained lower bound 
on the number of subharmonics of period kT, based on symmetry considerations 
and on pinching estimates, as in Sect. 5.2 of this article. 
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1 Fixed-Period Problems: The Sublinear Case 

With this chapter, the preliminaries are over, and we begin the search for periodic 
solutions to Hamiltonian systems. All this will be done in the convex case; that 
is, we shall study the boundary-value problem 



with H{t, ■) a convex function of x, going to -|-oo when — > oo. 
1.1 Autonomous Systems 

In this section, we will consider the case when the Hamiltonian H(x) is au- 
tonomous. For the sake of simplicity, we shall also assume that it is C^. 

We shall first consider the question of nontriviality, within the general frame- 
work of (Aoo, i?oo)-subquadratic Hamiltonians. In the second subsection, we shall 
look into the special case when H is (0, 6oo)-subquadratic, and we shall try to 
derive additional information. 

The General Case: Nontriviality. We assume that H is (Aoo, i?(x>)-sub- 
quadratic at infinity, for some constant symmetric matrices A^o and B^o, with 
Boo — ^oo positive definite. Set: 



X = JH'{t,x) 
x{Q) = x{T) 




7 : = smallest eigenvalue of B^ — A, 



■oo 



(1) 



A : = largest negative eigenvalue of J— + A, 



dt 



■OO 



(2) 



Theorem 21 tells us that if A + 7 < 0, the boundary- value problem: 



X = JH'{x) 
xifi) = x{T) 



(3) 



has at least one solution x, which is found by minimizing the dual action func- 
tional: 

^1 



-{A-'u,u)+N*{-u) 



dt 



(4) 



on the range of A, which is a subspace R{A)j^ with finite codimension. Here 



N{x) ■.= H{x)--{A^x,x) 



is a convex function, and 



(5) 



N{x) <-{{B^-A^)x,x)+c Mx 
Proposition 1. Assume H'{0) = and H{0) = 0. Set: 

S :=limM2N{x)\\x\\~'^ . 

a;— >0 

If J < —A < 6, the solution u is non-zero: 

x(t) ^0 yt . 



(6) 



(7) 



(8) 



Proof. Condition (7) means that, for every S' > S, there is some e > such that 

(9) 



\\x\\<e^N{x)<^-\\xf 



It is an exercise in convex analysis, into which we shall not go, to show that 
this implies that there is an r/ > such that 



f\\x\\<r]^N*iy)<^\\y\f 



(10) 



Fig. 1. This is the caption of the figure displaying a white eagle and a white horse on 
a snow field 

Since Ml is a smooth function, wc will have ll/iwill;,^ < rj for h small enough, 
and inequality (10) will hold, yielding thereby: 

Hhui)<~\\u,\\l + ^^\\u,f . (11) 
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If we choose 5' close enough to 5, the quantity {\ + jr) will be negative, and 
we end up with 

il){hui) < for /i ^ small . (12) 

On the other hand, we check directly that V'(O) — 0. This shows that cannot 
be a minimizer of ■0, not even a local one. So m ^ and u 7^ yl~^(0) = 0. □ 

Corollary 1. Assume H is and {aoo,boo)-subquadrafAc at infinity. Let ^1, 
...,^jv be the equilibria, that is, the solutions of H'{^) = 0. Denote by ojk the 
smallest eigenvalue of H" {^k), and set: 



:= Min {wi, . . . , Wfe} . 



If: 



2n 



boo < -E 



T 



T 



(13) 
(14) 



then minimization ofip yields a non-constant T -periodic solution x. 



We recall once more that by the integer part E[a\ of a G IR, we mean the 
a e ^ such that a < a < a -\- 1. For instance, if we take a^o = 0, Corollary 2 
tells us that x exists and is non-constant provided that: 



or 



—boo < 1 < 
2n 2n 



, 27r 27r 
TG ( — ,— 
w boo 



(15) 



(16) 



Proof. The spectrum of ^ is -\- aoo- The largest negative eigenvalue A is 
given by ^fc^ + Goo, where 



27r, 27r,, 

— fco + < < — (fco + 1) + Co 



Hence: 



ko = E 



T 



The condition 7 < —A < 6 now becomes: 



(17) 
(18) 



27r, 

doo < - jT^o - Cloo < W - 



(19) 



which is precisely condition (14). □ 

Lemma 1. Assume that H IS on ]R^"\{0} and that H"{x) is non-degenerate 
for any x ^ 0. Then any local minimizer xofip has minimal period T. 
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Proof. We know that x, or x + ^ for some constant ^ G IR^", is a T-periodic 
solution of the Hamiltonian system: 

X = JH'ix) . (20) 

There is no loss of generality in taking ^ = 0. So il){x) > ipi^) for ^ in 
some neighbourhood of x in W'^''^ [JR/TZ; IR^") . 

But this index is precisely the index irix) of the T-periodic solution x over 
the interval (0,T), as defined in Sect. 2.6. So 

irix) = . (21) 

Now if X has a lower period, T/k say, we would have, by Corollary 31: 

irix) = ikT/k{x) > kiT/kix) + k-l>k-l>l . (22) 

This would contradict (21), and thus cannot happen. □ 

Notes and Comments. The results in this section arc a refined version of 1980; 
the minimality result of Proposition 14 was the first of its kind. 

To understand the nontriviality conditions, such as the one in formula (16), 
one may think of a one-parameter family xt, T G (2nuj^^ ,2nb^^ of periodic 
solutions, xt{0) = xt{T), with xt going away to infinity when T — > 27rw~^, 
which is the period of the linearized system at 0. 



Table 1. This is the example table taken out of The TgXbook, p. 246 



Year 


World population 


8000 B.C. 


5,000,000 


50 A.D. 


200,000,000 


1650 A.D. 


500,000,000 


1945 A.D. 


2,300,000,000 


1980 A.D. 


4,400,000,000 



Theorem 1 (Ghoussoub-Preiss). Assume H{t,x) is {0,e)-subquadratic at 

infinity for all s > 0, and T-periodic in t 

H{t,-) is convex Vt (23) 

H{-,x) is T-periodic \/x (24) 
H{t,x) > n{\\x\\) with n(s)s~"^ — >■ oo as s — >■ oo (25) 

Ve>0, 3c : H{t,x) <^\\xf + c . (26) 
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Assume also that H is C^, and H"{t, x) is positive definite everywhere. Then 
there is a sequence Xk, fc e IN, of kT -periodic solutions of the system 

x = JH'{t,x) (27) 

such that, for every A; e IN, there is some Po € IN with: 

p > Po ^ Xpk ^ Xk . (28) 

□ 

Example 1 (External forcing^. Consider the system: 

X = JH'{x) + fit) (29) 

where the Hamihoniari H is (0, 6oo)-subquadratic, and the forcing term is a 
distribution on the circle: 

f=j^F + fo with F {M/T2Z; M^") , (30) 

where fo := f{t)dt. For instance, 

f{t) =Y.h^, (31) 

where 5k is the Dirac mass ai t = k and ^ G IR^" is a constant, fits the pre- 
scription. This means that the system x = JH'{x) is being excited by a series 
of identical shocks at interval T. 

Definition 1. Let Aac{t) and Boo{t) be symmetric operators mlR^", depending 

continuously on t £ [0,T], such that Aoo{t) < -Boo(^) for all t. 

A Borelian function H : [0,r] x H^" — > IR zs called {Aoo, Boa)-subquadratic 



at infinity if there exists a function N{t, x) such that: 

H{t,x) = ^{A^{t)x,x)+N{t,x) (32) 

Vt , N{t, x) is convex with respect to x (33) 

N{t,x) > n{\\x\\) with n{s) s~^ ^ +oo as s +oo (34) 

3cGlR : H{t,x) <^{B^{t)x,x)+c Vx . (35) 



If A^{t) = ttool and Boo{t) = bool, with aoo < 6oo G IR; we shall say that 

H is {aoo,boo)-subquadratic at infinity. As an example, the function \\x\\ , with 
1 < OL < 2, is {0, e)-subquadratic at infinity for every e > 0. Similarly, the 
Hamiltonian 

H{t,x) = ^k\\kf + \\xf (36) 
is {k, k + e)-subquadratic for every e > 0. Note that, if k <0, it is not convex. 
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Notes and Comments. The first results on subharmonics were obtained by Rabi- 
nowitz in 1985, who showed the existence of infinitely many subharmonics both 
in the subquadratic and superquadratic case, with suitable growth conditions 
on H'. Again the duality approach enabled Clarke and Ekeland in 1981 to treat 
the same problem in the convex-subquadratic case, with growth conditions on 
H only. 

Recently, Michalek and Tarantello (see Michalek, R., Tarantello, G. 1982 and 
Tarantello, G. 1983) have obtained lower bound on the number of subharmonics 
of period kT, based on symmetry considerations and on pinching estimates, as 
in Sect. 5.2 of this article. 
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